[1] "C"

Consideraciones iniciales

  • Variables independientes consideradas: Final.Glucose, NEFA, Cholesterol, Triglycerides, Glycerol, Lactate, Insulin, Glucagon

  • Variables respuesta consideradas:
    • De cada órgano: rna (mRNA), prot (proteína) y act (actividad)
      • rnaliver, protliver, actliver
      • rnakidney, protkidney, actkidney
      • rnamuscle, protmuscle, actmuscle
  • Los tamaños muestrales actuales son:

Group rnaliver protliver actliver rnakidney protkidney actkidney rnamuscle protmuscle actmuscle
Fasting.HFD 10 10 9 5 10 5 5 7 9
RF30min.HFD 11 11 9 5 11 5 5 7 9
RF3h.HFD 10 10 8 4 10 4 4 6 8
Fasting.SD 9 9 9 5 9 5 9 7 9
RF30min.SD 9 9 9 5 9 5 9 7 9
RF3h.SD 11 11 9 5 11 5 11 7 9

Teniendo esto en cuenta, a continuación se presenta:

  • Boxplot para el peso y la ingesta calórica, para comprobar que no hay sesgos experimentales inesperados.

  • Gráficos 3D que enfrenta mRNA, proteína y actividad de IDE, para comprobar visualmente si hay algún tipo de patrón o separación evidente entre las condiciones experimentales.

    • Anova multivariante (MANOVA) y gráficos de interacción que comparan el efecto de las condiciones experimentales en el mRNA, proteínas y actividad.
  • Regresiones parciales que analizan el efecto de las 8 variables independientes sobre cada variable respuesta (mRNA, proteína, actividad), basada en Análisis de Componentes Principales (PLSR-PCA, partial least square regression - principal component analysis).

Group vs Calories and Weight

El peso parece distribuirse uniformemente y de forma esperable entre las condiciones experimentales. La ingesta calórica de los ratones HFD.RF30min tiene una dispersión elevada en comparación con el resto de grupos, quizá sería conveniente explicar el posible origen de esta variabilidad.

mRNA vs Protein vs Activity

Se representa la relación entre la cantidad de mRNA, proteína y la actividad de IDE, en cada uno de los tres órganos. Los gráficos pueden rotarse. Buscad si visualmente ocurren cosas que se esperan (si hay más proteína cuando hay más mRNA, o más actividad), o si hay alguna agrupación que diferencie claramente una condición experimental de otra. Esto se cuantifica mejor en el siguiente apartado “MANOVA”.

MANOVA

A continuación se presenta el análisis conjunto de las variables respuesta (mRNA, proteína y actividad) y las diferencias significativas detectadas entre condiciones experimentales.

Para ello se ha llevado a cabo un Anova Multivariante o MANOVA, que ha dado como resultado que hay diferencias en las variables debidas a las condiciones experimentales (p < 0.01). Podéis ver dónde están esas diferencias en los gráficos de interacción que aparecen debajo, que representan las medias y las desviaciones:

  • Columna 1: Hígado
  • Columna 2: Riñón
  • Columna 3: Músculo

  • Fila 1: mRNA
  • Fila 2: Proteína
  • Fila 3: Actividad

Comparad aquellos intervalos que no se solapen, son los responsables de las diferencias significativas. Las comparaciones hay que hacerlas en dos sentidos:

  • Dentro de cada estado nutricional, se compara HFD (rojo) vs SD (azul)</span)
  • Dentro de cada tipo de dieta, se comparan estados nutricionales (rojos vs rojos y azules vs azules).

Para explicar esto en el artículo, sugiero señalar únicamente el p-valor del test omnibus (el MANOVA general que he comentado antes) y mostrar los gráficos, sin p-valores adicionales de comparaciones múltiples, a menos que los pidan los revisores. A lo sumo, para cada gráfico, se podría dar un p-valor adicional que indique si al menos una de las medias es diferente a las demás.

PLSR-PCA

Este análisis permite analizar la influencia conjunta de las variables independientes fisiológicas sobre las variables respuesta, y cómo afectan a la organización de las condiciones experimentales.

No devuelve p-valores porque es un análisis gráfico, que permite representar los datos en nuevas dimensiones que maximizan la variabilidad capturada.

En los siguientes apartados os explico el procedimiento y cómo interpretarlo. El resultado final del análisis es el apartado “Plot PLSR”.

Scale data

Primero, se escalan todas las variables, es decir, se hace que todas tengan media = 0 y desviación = 1. Esto es imprescindible para que las diferentes unidades de medida de las variables no sesguen los resultados del análisis. Podéis ver en el boxplot que las ditribuciones resultan simétricas, y todas están centradas en 0.

Este apartado no tiene mayor interés para vosotros.

Plot profiles

En este gráfico se muestra lo que se va a analizar con la técnica PLSR-PCA. Para cada condición experimental, tenemos un perfil de variables independientes (a la izquierda de la barra vertical) y tres perfiles de variables respuesta, uno por cada órgano (derecha de la barra vertical). Cada línea une los valores de una misma observación (ratón).

La simple inspección visual ya permite apreciar algunas diferencias, por ejemplo, para la dieta estándar SD, en la condición de Fasting toma protagonismo la variable NEFA y disminuyen los valores de Glicerol, Lactato e Insulina, sin embargo, este patrón parece invertirse para las condiciones de RF30min y RF3h.

El objetivo es ser capaz de analizar cómo estos perfiles de variables independientes, respuesta y condiciones experimentales se interrelacionan unos con otros. Eso se hace en el siguiente apartado “Plot PLSR”.

Plot PLSR

Se presentan a continuación 9 biplot. Un biplot es una representación, en el mismo espacio dimensional, de la correlación entre:

  • Variable respuesta (rna, proteína, actividad) y variables fisiológicas. Su efecto se representa mediante flechas: la flecha morada para la variable respuesta, las verdes para las fisiológicas.
  • Las condiciones experimentales, mediante puntos que reflejan los datos reales y un código de colores para identificar el grupo al que pertenecen.

Cada columna corresponde a un órgano, y cada fila a una variable respuesta de IDE.

Un ejemplo de interpretación, para el primer biplot, sería:

  • El biplot muestra la correlación entre los niveles de mRNA y las variables fisiológicas para el hígado.
  • Las variables (flechas verdes) que más explican los niveles de mRNA (flecha morada) son:
    • El glucagón, correlación inversa intensa (se alinea con la dirección de la flecha pero en sentido contrario). Esto implica que niveles altos de mRNA de IDE en el hígado correlacionan con niveles bajos de Glucagón.
    • Los triglicéridos y el colesterol, correlación directa leve o moderada (se alinean débilmente con la dirección de la flecha y en el mismo sentido). Esto implica que niveles altos de mRNA de IDE en el hígado correlacionan con niveles altos de triglicéridos y colesterol, en algunos casos.
  • El resto de variables apenas correlacionan con los niveles de mRNA (dirección casi perpendicular de las flechas verdes respecto a la morada). En concreto, el colesterol parece no participar de ninguna forma (flecha muy pequeña, casi sin separarse del centro del gráfico)
  • En cuanto a las condiciones experimentales:
    • Se observan dos cluster, bien diferenciados en función de Fasting o no Fasting.
    • El estado de Fasting, para cualquier dieta, correlaciona poco con los niveles de RNA. La variable que más les afecta es NEFA.
    • Los estados de no Fasting (RF30min y RF3h, para cualquier dieta), parecen más definidos por las variables Final.Glucose, Insulin y Lactate. Sus niveles de mRNA son mayores y más variables que los de Fasting (los puntos se orientan en el sentido de la flecha morada, pero están muy dispersos)

También es importante detectar cambios secuenciales entre los biplot. Por ejemplo, en el caso del hígado, se observa cómo las variables Final.Glucose, Insulin y Lactate pasan de correlacionar muy poco con la cantidad de mRNA, para aproximarse cada vez más a la flecha morada en proteína y más aún en actividad. Es decir, en hígado, la actividad de IDE está sobre todo influenciada por estas tres variables, lo que no ocurre a nivel de mRNA y proteína. También se observa que se mantiene la separación entre Fasting-no Fasting, y que en general el resto de variables no modifican su influencia.

Esta misma interpretación habría que hacerla para cada órgano y biplot, centrándose especialmente en los patrones que resulten llamativos y que puedan explicarse a nivel fisiológico.